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The meson (pion and kaon) form factor is calculated in the perturbative framework with alter- 
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native forms for the running coupling constant and the gluon propagator in the infrared kinematic 
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region. These modified forms are employed to test the sensibility of the meson form factor to the 
nonperturbative contributions. Its is a powerful discriminating quantity and the results obtained 
with a particular choice of modified running coupling constant and gluon propagator have a good 
agreement with the available data, for both mesons, indicating the robustness of the method of 
calculation. Nevertheless, nonperturbative aspects may be included in the perturbative framework 
' of calculation of exclusive processes. 
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I. INTRODUCTION 



The study of the infrared (IR) limit of the constant coupling of the Quantum Chromody- 
namics (QCD) had attracted much attention in the last years, as well as the infrared finite 
form of the gluon propagator. The behavior of the running coupling constant claimed for 
several theoretical and experimental studies indicates a frozen value in this kinematical re- 
gion The question of the infrared form of Green functions of QCD is still a controversial 
field of research (for a review see 0, [|). Due to the different methods employed to analyze 
the Green functions, for example, the Dyson- Schwinger equations (DSE) and simulations in 
the lattice field theory (LFT) (as well as the approximations used to avoid several difficulties 
found in these methods) we have distinct different forms of Green functions. However, as 
we will show, the combination of the results from LFT and solutions of DSE restricts the 
form of the gluon propagator. 

In a recent paper [|l[], Aguilar, Mihara and Natale perform a study of the running cou- 
pling constant in the infrared region, using the dependence of the meson form factors on the 
coupling constant. To calculate the pion form factor, a mixed framework between pertur- 
bative and non-perturbative physics is employed with the following assumptions: freezing 
of the coupling constant and finite gluon propagator in the IR region, although using the 
perturbative scheme of calculation of the meson form factor ||. 

Nevertheless, the use of a frozen running coupling constant to describe QCD exclusive 
processes, particularly the meson form factor, is not a novelty. Ji and Amiri || and Brodsky 
et al. || use a modified coupling constant to study the form factor and the related reaction 
77 — ► 7r + 7r~. This modified coupling constant is based on the work by Cornwall @ in which 
a gauge independent set of Dyson- Schwinger equations is solved, giving a coupling constant 
finite in the infrared and a gluon propagator with a running mass. 

However, there are another approaches to include infrared contributions. In ||, Maris 
and collaborators perform a full non-perturbative calculation of the meson form factor. The 
amplitude for the vertex Trrcy (or KKj), which is related with the form factor, is obtained 
using the following ingredients: the Bethe-Salpeter amplitudes for the scattering kernel qq 
(solutions of the Bethe-Salpeter equation); the solution of the DSE for the quark propagator 
in the rainbow truncation and an ansatz for the quark-photon vertex based in a Bethe- 
Salpeter equation in the ladder truncation. This approach gives a good description of the 



2 



experimental data for the low momentum region for both form factors (pion and kaon) as well 
as another observables (see last paper of ||). In another work ||, the authors employed a 
light-front Bethe-Salpeter model, where the quark, instead of the gluon as the model above, 
has a dynamical mass, giving a good agreement with the experimental data for the pion 
form factor. 

The meson form factor in turn, has many attempts of description in pure perturbative 
QCD UlQ| . For example, Stefanis UTTJ uses an unified factorization scheme that includes 



logarithmic corrections (which origin is the gluon emission) and power correction whose 



origin is non-perturbative. Melic et al. [0 calculate the pion form factor in next-to-leading 



order (NLO) in perturbation theory Yeh [HJ also made a NLO calculation and found a 
fittable expression for the pion form factor which includes non-perturbative effects. However, 
these attempts have only validity where the perturbation theory is valid, in other words, 
when the momentum transfer is large. 

In this article, we extend the results of |l| for the pion form factor for a different set 
of coupling constants and propagators and compare them with more recent experimental 
data. Also the robustness of the model is tested with the kaon, a more massive meson than 
the pion. As we will show, the model employed describes both, pion and kaon form factors 
with good agreement with the available data for the kinematical region of exchange of low 
momenta. The article is organized as follows: in the section |TJ, we will review the model 
used to calculate the meson form factor. Next, in section we present the modified meson 
form factor calculation scheme with the main features of the non-perturbative gluons and 
the frozen running coupling constant and their relation. Finally, we present the results and 
conclusions of the work. 

II. THE MESON FORM FACTOR CALCULATION 

The meson form factor is given by the factorized expression in perturbative QCD (pQCD) 
form [g, [TJ, [Tj| (see figure [I]) 

F M (Q 2 ) = dx dycf)* M (y,Q y )T H (x,y,Q 2 )(t)M(x,Q x ), (1) 
Jo Jo 

where Q z = min(,2, 1 — z)Q (z = x, y) and Q is the 4-momentum transfer by the photon. The 
equation above can be seen as a convolution of the initial and final states, represented by the 
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quark amplitude distributions 4>m(z, Qz) (obtained in a non-perturbative calculation) with a 
hard (perturbative) scattering amplitude, T H (x,y,Q 2 ). The quark amplitude distribution is 
interpreted as the amplitude to find the quark or antiquark within the meson with fractional 
momentum z or 1 — z, respectively. 




T H (x, V, Q) 



FIG. 1: The LO Feynman diagram for the meson form factor, where 4>m(x, Q 2 ) is the meson wave 
function. 



The function T H (x, y, Q 2 ) was calculated in leading order (LO) in |I6[ and in next leading 
order (NLO) in [JT3]. The LO expression is given by (fig. [I]): 

. 2 , 64tt (2a s {(l-x)(l-y)Q 2 ] 
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where a s (Q 2 ) is the running coupling constant. 

The quark amplitude distributions (Pm(x,Q 2 ) are usually calculated in a certain energy 
scale (Qo) by QCD sum rules |T7| and then evaluated in an arbitrary energy scale solving a 



Bethe-Salpeter-type evolution equation [|, p~q] , which solution is 
Mx, Q 2 ) = x(l -x)J2 Ci m (2x - 1) 

n=0 L««IV0JJ 

where Ci 3/,2 ^(2x — 1) are the Gegenbauer orthogonal polynomials |TI| 
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where A^ s is the non-singlet anomalous dimension given by 
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and /?o = 11 — | n / is the number of flavors). To obtain the factors <f>n(Qo), given by 
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(6) 



we use the quark amplitude distributions in a given scale of energy. From [H], we have at 
Qo = 500 MeV for the pion and the kaon, respectively, 

fa(x,Q 2 Q ) = ^ix{l - x)(2x - l) 2 , (7a) 
<Pk(x, Qq) = [0.6(2a; - l) 2 + 0.25(2a; - l) 3 + 0.08] , (7b) 
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where the normalization factor Jm is given by 



dx <j) M (x,Ql) = -0=, (8) 



and thereby f n = 93 MeV, f K = 112 MeV. 

Using the orthogonality relation of the Gegenbauer polynomials, it is easy to find that 
the only non zero 0n (Qo)' s f° r the pion are 

_ c f n j(!r)/ ( o2\ a fir 



2^3' " 2 V " U/ 2^3' 
and for the kaon 

#'(QS) = 6^, (9a) 



(Qo) = , CTO = ^ , (9b) 



^ 2 v " u/ 5 2V3 ' V " U/ 7 2v /3 

therefore the exponents d n are 



32 

4 = 0, 4 = 5^, (10a) 

d 2 = 50 , d 3 = — — • (10b) 

99 - 671/ ' 495 - 30n/ v ; 



Finally, the quark amplitude distribution is, for the pion, 



(,,g») = ^(l-,)J, + [30(2,-l)"-6](|g)*|, (11) 



and for the kaon, 
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However, Dziembowski and Mankiewicz [19] use the constituent-quark model to calculate 
the quark amplitude distribution for the pion, obtaining 

m 2 1 f (xM + m) [(1 -x)M + m] 



i (DM) , 



X 



N exp 
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(13) 



where N is determined by the normalization condition Eq.(H)(iV = 0.622), M is the spin- 
averaged meson mass (M = 614.4 MeV for the pion), m is the constituent quark mass 
(m = 330 MeV) and (3 is a Gaussian parameter chosen as 460 MeV. We will test the changes 
in the pion form factor when this quark distribution is employed. 



III. THE MODIFIED MESON FORM FACTOR 

In |l[], the modifications introduced in the pion form factor by the change of the usual 
running coupling amplitude by the frozen infrared one have been tested, as well as the case 
of the change of the gluon propagator. In this case, the hard scattering amplitude can be 
written as 

647r f 2 - - 1 1 
T H {x,y,Q 2 ) = — |-a.(* a )Z>(* a ) + -& s (f)D(f) j , (14) 

where k = (1 — x)(l — y)Q 2 , p 2 = xyQ 2 and D(Q 2 ) is the gluon propagator, and a s is the 
modified running coupling constant which expression is given below. 

Therefore, the expression for the meson form factor from Eq.(|l]) and (|I4]) is 

F M {Q 2 ) = dx dy4)* M {y,Qy)f H (x,y,Q 2 )4>M(x,Q x ), (15) 
Jo Jo 

where <pM is the meson wave function with the modified running coupling constant. 

The infrared form of the gluon propagator is still a controversial aspect (for a review for 
solutions and methods see and references therein). The general formula for the gluon 
propagator is given by (in the Landau gauge): 

K^)= 5ab (^- q -f-) D (^ ( 16 ) 

where the usual perturbative propagator is given by D(q 2 ) = 1/q 2 which diverges when 
q 2 — > 0. In the literature, there are several different forms for the gluon propagator, with 
different behaviors for the infrared region: finite, zero and more divergent that 1/k 2 , each 
one with its advantages and inconveniences. The reasons of these different behaviors are 
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the methods and approximations employed to obtain the gluon propagator. The more 
popular methods are the Schwinger-Dyson equations and the simulations from lattice field 



theory. The most recent results in the last method |20| , |2l|j discard the solution of the type 
l/k 4 and therefore we will not use it in this work. In spite of the controversy, the non- 
perturbative gluon propagator was used successfully in many phenomenological applications 



as, for example, proton-proton scattering and elastic production of vector mesons 




In this work, we use the following gluon propagators: 

• calculated by Cornwall @ using the gauge independent pinch technique, given by 

D -(f > = T+iiijSSj. (17) 

where M 2 (q 2 ) is a dynamical gluon mass term, 

-12/11 

(18) 

l QCD / / 

with rrig = 500 ± 200 MeV. With this propagator, Q obtains the best description for 
the pion form factor. Its features include the correct ultraviolet behavior (according 
to the renormalization group) and the dynamically generated mass M 2 (q 2 ). 

from Habel et al. ||24| , calculated using an approach that employs the same features of 
the perturbative theory, which implies a simplified set of Schwinger-Dyson equations, 
which solution is given by 

Dn(q 2 ) = — ^, (19) 

We should note that this propagator goes to zero when q 2 — > 0, unlike the other 
propagators employed in this work, which have a non-zero value in this region. 



Alkofer et al. J25| solve a set of Dyson- Schwinger equations for the QCD Green's 



functions and the gluon propagator found by this manner can be fitted M by 

DAi, 2 ) = (20) 

where b = 3.707 and c = 0.603. 
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Atkinson and Bloch |26| solve a set of Dyson-Schwinger equations for the gluon and 
ghost propagator, improving the solution found by |2IJ. The asymptotic infrared 
behavior of the gluon propagator is 



if) - j 2 |^ (g 2 ) 2K f^ 2 )^ } . (2i) 



where A = I, k = 0.769475, a x = -10.27685, p = 1.96964, = 1, / 2 = 0.408732, 
/ 3 = -0.761655. 



from Gorbar and Natale p7j|, that use the operator product expansion (OPE) to 



relate the gluon and quark propagators with their respective condensates, the gluon 
propagator obtained is 

^G N (g 2 ) = — -i , (22) 



where pr = 0.61149 GeV is a parameter fixed by the gluon condensate |27] and 



X =0.9666797. We should note that this propagator interpolates two different propa- 
gators: the pure massive and one similar to the Habel solution. 

The modification of the gluon propagator is closely related with the modification of the 
running coupling constant in the infrared region. As calculated by Cornwall , the coupling 
constant is frozen in the low momentum by the addition of the massive term, Eq.flT8|), 



9 ' 

~f? 



where Aq CD is the QCD scale parameter and (3q — 11 — ^n/ is the first coefficient of the beta 
function and n/ is the number of flavors. In the case of the other propagators, we consider 
the terms in the denominator as massive terms and we will substitute it in the Eq.(|23[) as 
well as in the Cornwall propagator, as shown below 



The Habel propagator gets 



a W(a 2 ) = — (24) 



,l QCD 

with 



8 



The Gorbar-Natale propagator gets 

„(GN) (g2) 



All 



a) 



A) In 



with 
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(25) 



M^ N (g 2 )=/i 2 g ^-g 2 ) + ^f% 
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However, in the case of the propagators given by the equations (£(]) and fl2T|) , the cou- 
pling constant comes from the solutions of the Dyson-Schwinger equations for the Green's 
functions. For the case of the solution of the Alkofer et al, the coupling constant can be 
fitted by the following formula M 



"I 

Oil 



q 2 < 0.31 GeV 2 

0.31 < q 2 < 1.3 GeV 2 

q 2 > 1.3 GeV 2 



(26) 



with 



0.2161 + 9.2621 exp 



a h s (q 2 ) = 1.4741 + 8.6072 exp 



2 (g 2 - 0.0297)^ 
' (0.6846) 2 
g 2 - 0.1626\ 



0.3197 



a c s (q 2 " 



1.4978 



(27) 



ln(1.8488g 2 ) 

The solution of Atkinson and Bloch |26[ is given analytically in the asymptotic regions 



k=l 



2 \ k P 



n 2 J 



(28) 



where v = 0.912771, ft 2 = 0.1864754, b x = -1, b 2 = 0.760753, 6 3 = -0.370785, and 



Alt 



41og( 



(29) 



A QCD' 



where A^ CD = 0.06802 GeV 2 . 

The procedure is the following: using the Eq . (|15D , we substitute in the hard scattering 
amplitude Tjj the different gluon propagators in the order above as well as the running 
coupling constant, a s , to verify the changes in the form factor. A comparison with the 



perturbative fit of Yeh |13| , as well as, with the full non-perturbative calculation of Maris 
and collaborators M is also performed. 
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FIG. 2: The pion form factor calculated using the modified running coupling constant and gluon 
propagator from different works, indicated in the table above. The number of flavors is zero. The 
labels CZ and DM refer to the eqs.(ll) and (13) for the quark amplitude distributions. The solution 



of Maris et al. from [g]. The Yeh fit from eq. (31). Data from [g8|, pSj ]. 
IV. RESULTS 

The results obtained for the pion form factor in the flavorless case (rif = 0) are displayed 
in the figure |2|. A remarkable feature of the result is the significant difference between 
the distinct coupling constants and their respective gluon propagators. The origin of this 
difference is the strong dependence of the form factor in the running coupling constant. The 
results are also consistent with the infrared behavior of the coupling constant, for example, 
in the case of the Habel propagator (Eq. (^)) when Q 2 —>■ the massive term diverges, 
which implies a divergent behavior. With the other propagators, the behavior is distinct 
because the coupling constant is finite when Q 2 — > 0, giving a finite result. 

As already pointed out in Q, the best description for the data is given by the Cornwall 
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FIG. 3: The pion form factor calculated as in figure |2| with nj = 3. 

propagator with m g = 300 MeV while the other choices for the propagator listed above do 
not give a reasonable result in comparison with the available experimental data. 

When a more realistic number of flavors, n/ = 3, is used, the result for the form factor 
does not present a significant change, as expected in U and displayed in fig. (|3|), however 
the result with three flavors is better than with no flavors. 

As the number of flavors, the change in the quark amplitude distribution, Eq. fllTD and 



(13|) does not make a significant modification in the form factor. 

We display in the figures (|2]), (|3]) and (f|), the full non-perturbative solution for the pion 
form factor obtained by Maris and collaborators ||, only in the low momentum transfer 
region {Q 2 < 0.8 GeV 2 ). As pointed out in the Introduction, this solution gives a very good 
description of the low momentum data. Our best result, using the Cornwall's propagator 
and m g = 0.3 GeV, is below of this prediction although with the same global behavior. If 
we diminish the smaller massive parameter m g , we can obtain a better fit to the data (for 
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FIG. 4: Plot of Q 2 F 7T versus Q 2 using the same coupling constant and propagators used in the fig. 
(|2|), with quark amplitude distribution from Eq. (|i~l| ) and rif = 3. 

both mesons), but this value is outside of the range of the values for this parameter, which 
provides good description of another processes [p2j ^3 . 



We also display a fit with the experimental points, based on the approach of |13[ . In that 
paper, the author uses the collinear expansion to analyze the NLO power corrections to the 
pion form factor, and founds the following expression for the pion form factor 



FM) - ma ^ f ' (l - ) ■ (30) 



where Q 2 S is the effective energy scale and J{Q 2 eS ) is called a jet function, introduced to 
absorb the infrared divergences from the NLO corrections. In the same paper, it was made 
a phenomenological fit to the data for the pion form factor with the following form 



*?«A = £(i-!). (3D 
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FIG. 5: The kaon form factor calculated as the pion form factor, but using only the quark amplitude 
distribution from Eq. (|T^). The phenomenological fit is given by Eq. (^). Data from |?l]] 



In [|T^J, the best values for A and B to fit the above formula to the experimental data from 



3TJ are A = 0.46895 and B = 0.3009. 
The phenomenological fit has a divergence when the momentum transferred by the photon 
goes to zero, in opposition to the results that employ the non-perturbative propagators and 
a finite infrared running coupling constant, and therefore does not give a good description 
for the data in low Q 2 . 

In order to compare with the previous results found in the literature, we plot Q 2 F n (Q 2 ) 
as a function of Q 2 , as shown in the fig. [|. In 0, its prediction does not go to zero, as 
in our result, neither describes the data in the region of low Q 2 . In p9 |, a comparison of 
predictions found in the literature for the pion form factor is made and the best result of 
this study (Cornwall with m g = 300 MeV) still is a good result. 

In the case of the kaon, we used the same procedure, although using only the quark 
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FIG. 6: The kaon form factor calculated as in figure (||), with 



n, 



amplitude distribution of the Eq. fll^) in the Eq.fll4l). The result of the flavorless case (n/ — 0) 
is displayed in the fig.® while the case rif = 3 is shown in the fig.(||). As shown in the 
figures, the global behavior is the same as in the pion form factor, in other words, the best 
description for the data is given again by the Cornwall's choices for the coupling constant 
and gluon propagator with m g = 300 MeV. 

In the figures @ and (|j) we present a result, based in |13|], with the same functional form 



Fk\Q 2 ) = 7V2 i-nj). ( 32 ) 



of the pion, but now applied to the available data for the kaon, 

Q 2 \ Q 2 

since the only dependence of the meson in the Eq.(j30"l) is the normalization factor J'm- 
We compare Eq. (|32l) to the data available for the kaon and found for the coefficients 
C = 0.0594695 and D = 0.0130963 with X 2 = 1-93314. 

As the data available for the kaon restricts to the region of low Q 2 , the predictions of 
different models cannot be tested in the kinematic region of high Q 2 for the moment, in 
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FIG. 7: Plot of Q 2 Fk versus Q 2 using same coupling constant and propagators used in the fig. 
(|5|), with quark amplitude distribution from Eq. ( |T^ ) and n/ = 3. 

opposition to the case of the pion. 

For the pion, as well as for the kaon, when the coupling constant and the propagator 



of Atkinson and Bloch [26[] work is employed, due to the absence of a full analytical form 
for the coupling constant and for the gluon propagator, we use only the infrared analytical 
forms, exception to the case of the fixed scale of energy (see Eq.([7a|)), when we use the ultra- 
violet analytical result. There is a difference in the full numerical result and the analytical 



asymptotic expression for the coupling constant (see the first paper of |26[]), but it is not 
significant for our result. 
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V. CONCLUSIONS 



In this work, we calculated the pion and kaon form factor following the approach of 
j ij obtaining a reasonable agreement with the data available (for the pion [^8|, ^ |3(J and 
kaon [|3T]) in the case of the Cornwall's gluon propagator |7|] (with parameter m g = 300 GeV) 
and the running coupling constant that is frozen in the small momentum transfer region, 
avoiding infrared divergences. We point out that a smaller mass gives a better result, but 
its value is out of the mass interval found in previous results for another processes P5[ . 
The other propagators (Cornwall with m g = 700 GeV, Habel el al. [E3J], Alkofer el al. [25 



Atkinson and Bloch and Gorbar and Natale |27|[ ) give results in disagreement with the 
experimental results. The significant difference between the results has its origin in the well- 
known meson form factor sensibility to the running coupling constant. Since we modified 
the IR behavior of the a s , this difference is expected as in [jlj, where the same behavior is 
observed. The case with a non-zero number of flavors, rif = 3, does not have a significant 
difference in the results. The exchange of the quark amplitude distribution in the pion case 
also does not modify significantly the form factor according to ||. The model employed 
also describes the available experimental data in the low Q 2 region for both (pion, kaon) 
mesons form factors, in opposition to the results || that employ only a frozen form of the 
running coupling constant, without the finite infrared form for the gluon propagator, but 
agrees with the results shown in |29[ for the pion form factor. In [JT2, 5^], restrictions are 
made concerning the applicability of the frozen coupling constant and the modified gluon 



propagator. In |32j the pion form factor is calculated using the same ideas of this work, 
although employing a different pion wave function and a gluon propagator with fixed mass, 
resulting for Q 2 F 7T (Q 2 ) a smaller value than the experimental one by a factor of 10. In 
comparison with the pure non-perturbative calculation of |§ and the NLO perturbative 



calculation of | 13fl , our results for the pion (except by a normalization factor, given by 
adjusting the massive parameter) is an interpolation of the above results, having the fit 
characteristics of both calcutations. The good agreement with the experimental data of this 
work indicates that it is possible to include non-perturbative effects in exclusive processes 
in QCD without dramatic changes in the perturbative scheme of calculation of observables. 
Otherwise, the inclusion of non-perturbative effects is still a field in which much research is 
required. 
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